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Abstract 

Kinetics of magnetic flux in a thin mesoscopic ring biased by a strong external mag- 
netic field is described equivalently by dynamics of a Brownian particle in a tilted wash- 
board potential. The 'flux velocity', i.e. the averaged time derivative of the total magnetic 
flux in the ring, is a candidate for a novel characteristics of mesoscopic rings. Its global 
properties reflect the possibility of accommodating persistent currents in the ring. 

1 Mesoscopic rings: two-fluid model 

Mesoscopic devices have attracted much theoretical and practical attention because they are 
promising for implementation in ultra- small hybrid elements to test quantum information the- 
ory [1]. A large class of such devices is based on ring structures, i.e. the Aharonov-Bohm 
topology. Such a class contains both superconducting (SQUIDs) and non- superconducting de- 
vices. 

In this paper we study selected kinetic aspects of persistent currents which can be observed 
in normal metal, semiconducting rings or cylinders and, as probably the most famous exam- 
ples, in carbon nanotubes or nanotori. We focus our attention on kinetics of magnetic flux in 
the presence of a strong external static magnetic field. We show that it can be modeled in the 
same way as the dynamics of a Brownian particle moving in a biased washboard potential. 
Here, the analog of the position of the Brownian particle is a total magnetic flux. We show 
that the time derivative of the magnetic flux, i.e. the flux velocity (if we recall the analogy to 
the dynamics of the Brownian particle) depends strongly on the ability of accommodation of 
persistent currents by the ring. 

Persistent currents are equilibrium currents flowing in the Aharonov-Bohm systems which 
are small enough to preserve phase coherence of electrons [0 [3J . In ideal samples at the van- 
ishing temperature T = 0, all electrons are the carriers of such a current. It is not the case 
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at non-zero temperatures T > 0, when some of the electrons are no longer coherent and are 
a source of the 'normal' Ohmic current. Let us consider now a mesoscopic ring placed in a 
uniform magnetic field B in the 3-dimensional space. Because of the self-inductance L, the 
electric current / will induce a magnetic flux in the ring. Therefore, the flux and the current 
in the ring are coupled according to the expression 

4> = <p e + LI = <p e + L[I coh + I dis }. (1) 

The flux <p e is induced by the external magnetic field B. The total current / is a sum of the 
coherent current I coh and the Ohmic dissipative current I dis . The coherent current is assumed 
to be a linear combination of the paramagnetic and diamagnetic contributions. This is related 
to occurrence, with a probability p, of the so called current channel with an even number of 
coherent electrons or an odd number of coherent electrons, with a probability 1 — p. Hence, 
with O = h/e, it reads 



T) = I J2M T /T*) jpsin (2™^ + (1 -p)sin 2vm (J- + ~ 
The amplitudes take the form 



'co/i *coh 

n=l 



(2) 



A n (T/T*) = 4 J eXp( , Wr/ ri, cos(nM*) • (3) 
v 1 ' tiT*1- exp(-2nT/T*) v ' 

The characteristic temperature T* is determined from the relation k B T* — A F / 2n 2 , where A F 
marks the energy gap, k F is the momentum at the Fermi surface and l x is the circumference of 
the ring. The parameter I is the maximal value of the persistent current at temperature T = 0. 
The dissipative current Idis is determined by the Ohm's law and Lenz's rule flU, 



his = i dis (<P,T) = ---g + r(t) , (4) 

where R is resistance of the ring, ks denotes the Boltzmann constant and T(t) describes ther- 
mal, Johnson-Nyquist fluctuations of the Ohmic current. This thermal noise is modeled by 
the Gaussian white noise of zero average, i.e., (T(t)) = and 5-auto-correlation function 
(r(i)r(s)) = 5(t — s). The noise intensity D = ~j R is chosen in accordance with the 

classical fluctuation-dissipation theorem. 

From equations (l)-(4), we get the Langevin equation governing the dynamics of the mag- 
netic flux |[6l : 

^ = ~ 0e) + W0, T) + sj^fm- (5) 

It can be rewritten in the form 



Rdt d(h V R 
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where the generalized potential W{4>) reads 
1 



w{4>) 



2L 



(0 2 - 20 e 0) 



^AJT/T\ , 
+ (polo — 7tz — \ P cos ( 27rn " 



n=l 



Inn 



p) COS 



2ixn 



+ 1/2 



•(7) 



Eq. © has been analyzed under various regimes [6J. In the following, we study specific regime 
of this system. 




Figure 1 : The rescaled generalized poten- 
tial ([8]) for p = 1/2 and temperature T = 
T/T* = 0.5. The potential approaches the 
washboard potential (flOl) when the exter- 
nal flux x e increases and the ratio x e / L\ = 
4> e / (I L) is fixed. This ratio is 1/5 and 
x e = 1, L x = 5 (solid line), x e — 3, L\ = 
15 (dashed line), x e = 5, L\ = 25 (dash- 
dotted line) and x e = 20, L x = 100 (dot- 
ted line). 




Figure 2: The washboard potential V(x) 
defined by Eq. (TTOb for p — 1/2 and tem- 
perature T = 0.5. For the rescaled cur- 
rent amplitude io > i c , the potential barri- 
ers appear. The critical value of the current 
amplitude i c is defined by the conditions 
V'(x) = and V"{x) = 0. For the pre- 
sented set of parameters i c ps 2. 



2 Flux-biased regime 

We intend to investigate the flux-biased regime which is defined in the following way 0: 
Let the external magnetic field B increases giving rise to increase of the magnetic flux e . 
Let additionally the self-inductance L increases. Formally, we perform the limit <p e — > oo 
and L — > oo in such a way that the ratio (j> e /L = I e is fixed. In this limit, the generalized 
potential W((f>) approaches a washboard form. Indeed, in Fig. 1, we present four forms of the 
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dimensionless generalized potential 

W 1 (x) ee = -L (a; 2 - 2x e x) 



WO 



— - {pcos (27rmr) + (1 - p) cos [2vm (x + 1/2)]} , (8) 
27m 

n=l 

with the dimensionless flux x = 4>/4>o, x e = 4> e /4>o an d the dimensionless inductance L\ = 
LI /(f) . We can notice that for the fixed ratio x e /L 1 = <fi e /(LI ) = I e /I = 1/5 and for 
x e = 20 and L\ = 100, the potential W\{x) is very well approximated by the biased washboard 
potential for large (but finite) number of periods of the coherent current I co h(<P, T). In the flux- 
biased regime, the Langevin equation (5) takes the form 

dx dV(x) 



+ V2D£{t), (9) 



dr dx 

where the dimensionless time r = t/r with r = (f) /RI e and the biased washboard potential 
V (x) (see Fig. 2) reads 

V(x) = V(x,T ) = —W 1 (x) 



" A (T ) 

x + ioY^ nK ° J {pcos(2vrnx) + (1 - p) cos \2irn(x + 1/2)1} (10) 

n=l 



27m 



with the rescaled current amplitude i = Io/I e . The rescaled zero-mean Gaussian white noise 
£(t) has the same 5-auto-correlation function as T(t). Its intensity D = k T , where the di- 
mensionless temperature T = T/T* and k = ksT* /(p I e is the ratio of thermal energy at the 
characteristic temperature T* to the energy of the flowing current I e induced by the elementary 
flux (f) . 

The Langevin equation © can be interpreted in terms of the overdamped motion of the 
Brownian particle in the washboard potential (flOl) . The periodic part of this potential is a 
'ratchet type' potential (8]|, i.e. V(x) does not posses the reflection symmetry. Let us notice 
that for p = 1/2, there is an additional periodicity V'(x + 1/2) = V'(x) presented in the 
system. 



3 Flux velocity 

In this paper we shall study the averaged (with respect to the noise realizations) stationary flux 
velocity (x) which is given by the formula flU: 

'V(l,T )-V(p,T 



N 



-i 



D 



-i 



(x) = N 
i 

dx 



1 — exp 

x+l 



dyexp 



D 

V(y,T )-V(x,T 
D 



(11) 
(12) 
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Figure 3: The averaged stationary flux ve- 
locity as a function of temperature T for 
several (both subcritical and supercritical) 
values of i and fixed p = 1/2. Inset: The 
corresponding washboard potentials (TTOl) 
at T = 1 are depicted for i = 2 (upper 
panel) and i = 10 (lower panel). 



Figure 4: The averaged stationary flux ve- 
locity as a function of T for three values 
of p and fixed i = 10. Inset: The cor- 
responding washboard potentials (flOl) at 
T = 1 (upper panel) and T = 2 (lower 
panel) are presented for three values of 
p = 0,0.25,0.5. 



The flux velocity is a function of the system parameters. The first parameter is the temper- 
ature. The second one is the current amplitude i , which reflects the ability of accommodating 
persistent currents. The third one, p E [0, 1], describes the structure of current channels as it 
is the probability of occurring current channel carrying even number of phase coherent elec- 
trons. Let us notice that, due to quantum size effects, persistent currents are always present in 
a sufficiently small system, i.e. there are always electrons maintaining their phase coherence 
when moving around the ring. The problem is if i is sufficiently large for those electrons to 
produce significant contribution to the total current flowing in the system at a given temper- 
ature. Numerical results show that upon inspection of the properties of the flux velocity (x) 
one can infer when the given ring is able to accommodate persistent current of a significant 
amplitude at a given temperature. In Fig. [3] we present the relation between the flux velocity 
and the temperature T for several different values of the current amplitude i . The general 
tendency is that in the presence of persistent currents the flux velocity is suppressed at the low 
temperature T or large %q. There are two classes of the systems split by the critical value i c of 
the current amplitude i , which determines the inflection points of the potential. This qualita- 
tive change is defined by a set of two equations: V (x) = and V"(x) = and is presented in 
Fig|2] For i > i c , barriers of the potential V(x) exist. For i < i c , the potential is a monotonic 
function of the flux x. Systems from the first class, with the supercritical amplitude i > i c , 
exhibit vanishing flux velocity for T — > 0. For the second class systems, with the subcritical 
amplitude io < i c , the flux velocity decreases but remains finite as temperature T — > 0. It is 
clear that in the formal limit i = 0, the flux velocity is constant, (x)u 0=0 = 1. The critical 
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value i c m 2 is estimated for the ring with p — 1/2. 

In Fig. 4, the probability p of a channel with an even number of coherent electrons is 
chosen as a parameter. For the system with statistically equal number of channels of both types 
(p = 1/2), the flux velocity is greater than for systems with one type of channels dominating 
over the other (e.g. p = 0). The results obtained for p = 1/4 and p = coincide with p = 3/4 
and p — 1 respectively and there is no way to distinguish with type of channels dominates in 
the ring. 

In order to quantify the effect, one can define the 'susceptibility', i.e. the temperature 
derivative of the averaged stationary flux velocity at fixed values of other parameters. Its mono- 
tonicity characterizes the possibility of obtaining persistent currents in the system. With this 
function one can associate a measure of the ability of accommodating persistent currents in the 
ring. This measure could be defined as a distance, in the sense of a metric in a function space, 
between the given (non-zero) susceptibility and the zero susceptibility. 

In conclusion, we have shown that performing suitable limiting procedure one can obtain 
new significant informations about persistent currents in mesorings. Investigations of global 
properties of the flux velocity can serve as an additional characteristics of mesoscopic rings. 
The perfect example is the monotonicity of the temperature derivative of the flux velocity 
or its asymptotic behavior at low temperatures which carries information about possibility of 
appearing persistent currents in the ring. 
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